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ABSTRACT 

Context. We come out in this paper with an astroclimatological study of meteorological data on relative humidity, dew point tem- 
perature, air temperature, wind speed, barometric air pressure, and sky cloudiness recorded at four Romanian locations (Baisoara, 
Rosia Montana, Semenic, Ceahlau) and Nordic Optical Telescope (NOT) located at Observatorio del Roque de Los Muchachos, in 
the Canary Islands. 

Aims. Long term trends of microclimates have been compared in order to identify the site-to-site variations. We have performed this 
analysis as part of a site testing campaign aimed at finding the best location for the establishment of a small Cherenkov telescope in 
Romania. The conditions at the Romanian sites have been compared to those of the Canary Islands considered as a reference. 
Methods. A statistical approach has been used for data analysis. Monthly and annual samples have been extracted from series of 
raw data for nighttime, daytime and entire day intervals. Percentage distributions of meteorological parameters, whose values exceed 
certain limits, have been derived. 

Results. Significant differences have been found between the Romanian sites and the NOT site. The comparison of the Romanian 
locations indicates Baisoara to be the best site for the establishment of the telescope, closely followed by Rosia Montana. As these two 
sites are both located in the Apuseni Mountains, we consider this area as the optimal place for performing astronomical observations 
in Romania. 
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1. Introduction 

The very high-energy (VHE) y-rays (100 GeV < E < 100 TeV), 
of cosmic origin, are best detected by the ground based Imaging 

' Atmospheric Cherenkov Telescopes (lACTs) (Buckley et al. 

; I2008i). M AGIC (Baixeras et al. i2003i), VERITAS (Weekes et 
aL l2002l l. H E.S.S . (Hinton et aL l2004l l. and CANGAROO-III 
(Kubo et al. 120041 ). the most important lACT facilities currently 

, in operation have reported scientific results with a huge im- 

■ pact on astrophysics, cosmology and particle physics. The future 

■ Cherenkov Telescope Array (CTA) (Hofmann & Martinez l2010l l 
with its next generation, highly automated lACT telescopes will 
open a new era of outstanding precision for gamma-ray astro- 

^ physics. 

All major lACT experiments perform observations at their 
sensitivity limit or in multi-wavelength campaigns for most of 
their operation time. If monitoring observations take place, the 
amount of time assigned is, by far, not sufficient. Under these cir- 
cumstances and considering the importance of the observational 
data, a global network of several small Cherenkov telescopes 
was proposed to be operated in a coordinated way for long-term 
monitoring of the brightest blazars - the Dedicated Worldwide 
AGN Reseai-ch FaciHty (DWARF) (Bretz et al. 2008). FACT 
(First G-APD Cherenkov Telescope) (Backes et al. i2011l) . the 
prototype telescope of this network has been prepared for opera- 
tions at the Roque de los Muchachos on the Canary Island of La 
Palma. 



To join the international efforts on understanding the physics 
of VHE y-rays, a small, ground based, imaging atmospheric 
Cherenkov telescope of last generation is intended to be installed 
and operated in Romania as a component of the DWARF net- 
work. 

A site assessment campaign has been initiated for the iden- 
tification of the best location where to establish the telescope. 
Technical, economical and social selection criteria were applied 
on a list of possible sites and four candidate locations have been 
selected (Baisoara, Rosia Montana, Semenic, Ceahlau) (Radu et 
al.|2010). An astroclimatological study based on a statistical data 
analysis of local meteorological parameters has been performed 
for these locations and the results are presented in this paper. The 
local microclimates at the Romanian candidate sites have been 
compared to that of Nordic Optical Telescope (NOT) operated 
on the island of La Palma (Canary Islands, Spain) (Djupvik & 
Andersen 2009 ). 

The international scientific community has sustained over 
the years large efforts to find the best sites for the operation of as- 
tronomical telescopes, but now it is for the first time that an anal- 
ysis with impact on the operation of a ground based Cherenkov 
telescope has been performed in Romania. 

The layout of the paper is as follows: in ^ we discuss 
the influence of the atmospheric conditions on the operation of 
Cherenkov telescopes; in Othe statistical characteristics of the 
data used for analysis are presented; in ^we include the astro- 
climatological results of this study. 
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2. Influence of the atmospheric conditions on the 
operation of Cherenkov telescopes 

As for most other types of astronomy, good atmospheric condi- 
tions are extremely important for the astronomical studies per- 
formed on very high energy gamma-rays with Cherenkov tele- 
scopes. The atmospheric conditions generally affect the quality 
of astronomical observations as well as the safety of the instru- 
ments and they are quantitatively evaluated by meteorological 
parameters. The observers monitor these parameters on a con- 
tinual basis in order to ensure a correct treatment of the astro- 
nomical data and the safe operation of telescopes (Cogan l2006l l. 

In the specific case of Cherenkov telescopes, it was pointed 
out that atmospheric changes on relative humidity, air tempera- 
ture, sky cloudiness and the influence of haze and calimaQ usu- 
ally not taken into account in simulations, may lead to a sys- 
tematic underestimation of the shower energy by ~ 15% (Mazin 
I2OO7I 1. 

As a function of the type of observations performed and the 
specific design of each telescope, the staff enforces safety limit 
values for the relevant meteorological parameters beyond which 
the telescope cannot be operated anymore. Poor weather is gen- 
erally considered as the major cause of the lost observing time. 
In comparison, the telescope failures and the engineering ac- 
count for a much smaller fraction of the lost time. This fraction 
was reported to be only ~ 3% in the case of VERITAS (Grube 
I2OO7I 1. 

In the following subsections we review the influence of the 
most significant meteorological parameters on the operation of 
Cherenkov telescopes. 

2. 1. Relative humidity and dew point temperature 

The relative humidity and the dew point temperature set the oc- 
currence of moist and water condensation on the coldest parts of 
astronomical instrumentation. If dew point temperature equals 
air temperature and if air cools, condensation appears. However, 
observing is dangerous for instrumentation even though conden- 
sation has not actually shown up, but there is a high risk that this 
phenomenon will occur This happens if the diff'erence between 
the air temperature and the dew point temperature is smaller than 
a variable upper limit (1 °C or 5 °C, Lombardi et al. i2009i) . 

As a result of condensation the light sensors and the elec- 
tronic equipment become wet and can be damaged. In case the 
condensation occurs on the optical surfaces which include the 
upper surface of the telescope reflector and the interior surfaces 
of the light concentrators, the image of the observed astronom- 
ical object degrades. In what concerns the relative humidity, a 
value of 90% was used by the NOT collaboration, by Lombardi 
et al. (2007 ), and by Murdin ( 11985D as a threshold beyond which 
the observations should be stopped and the telescopes turned off. 
Mazin (^007 ) and Jabiri et al. (I2000I I also considered that a rela- 
tive humidity level higher than 90% makes observing unfeasible. 

High levels of relative humidity also result in an increase of 
air conductivity. This is important for the light detectors located 
in the camera of a Cherenkov telescope. Especially the photo- 
multiplier tubes (PMTs) and their voltage supplies subjected to 
very high voltage during normal operation can experience elec- 
trical discharges which generate irreversible damages to the cir- 
cuits (Cogan|2006l gelik|2008l Theiling |2009] l. 



' fine dust or sand carried by winds from the Sahara Desert 



2.2. Wind speed 

Astronomical observations are not undertaken in windy weather 
due to dangers associated with the wind load on the optical re- 
flector of the telescope, on its support structure, and on the drive 
subsystem (Cogan 2006 ). In order to ensure that none of the sub- 
systems of the telescope are in danger, the telescope should be 
turned off and brought to the parking position if the wind speed 
exceeds an upper safety limit. This limit is telescope dependent. 
A survival wind speed limit is also defined for the parking posi- 
tion. 

The telescopes of VERITAS have positioner^ designed to 
operate safely at wind speeds up to 20 mph (8.9 m km 
/i ') and to survive in the parking position to wind speeds of 100 
mph (44.7 m .rVl61.0 km /i-') (Cogan , 2006. gelik2008). For 
MAGIC, it is considered that the subsystems of the telescopes 
are not in danger as long as the wind speed is lower than 40 km 

(11.1 m,ri) (Mazin 2007). 

Various other safety wind speed limits are discussed in the 
literature. The telescopes in operation at Observatorio del Roque 
de Los Muchachos (La Palma, Canary Islands), are closed if the 
wind speed exceeds 15 m 5"', whereas those at the site of Paranal 
Observatory (Atacama Desert, Chile) are closed for wind speeds 
in excess of 18 m (Lombardi et al. 2009 ). In a meteorological 
study performed for the Oukaimeden site in Morocco, 15 m i ' 
was also indicated as a typical maximum safe operating value 
for the wind speed (Murdin 1985; Jabiri et al. 2000). 

If high speed winds carry large amounts of dust the negative 
impact on the various hardware subsystems of a telescope is even 
worse. Besides that, dust is also a serious cause of decreasing 
transparency. The quantity of dust in the air depends on the alti- 
tude of the observing site, its proximity to a dust source, and the 
prevailing wind direction. From this point of view the Romanian 
sites located in areas covered with vegetation and away from 
dust sources are privileged. 

2.3. Air temperature 

The air temperature fluctuations impact on the instrumentation 
of Cherenkov telescopes and on the reconstruction of shower 
images. 

The PMT noise increases with temperature due increased 
thermionic emission whereas the PMT gain reduces with tem- 
perature. Cogan (2006) reported a gain decrease of (~ 1% / °C 
for these sensors. The temperature inside the camera of the tele- 
scope can increase as a result of the exterior temperature or be- 
cause of the heat dissipated by the electronic components lo- 
cated inside it (ex. preamplifiers). This is especially danger- 
ous while performing maintenance operations during daytime in 
warm seasons. As a safety precaution, exhaust fans and temper- 
ature sensors are usually installed inside the telescope camera 
(gelik|2008l Theihng|2009). 

On the other hand, temperatures below °C generate freez- 
ing. In order to avoid dew and snow effects, the mirror tiles of 
the large area reflector of MAGIC telescopes contain heating el- 
ements (White I2007] ). 

For Cherenkov telescopes, the main contribution to atmo- 
spheric attenuation comes from Rayleigh scattering, which de- 
pends on the chemical composition and the density of air 
molecules. The density of air molecules depends in its turn, on 
air temperature and barometric air pressure (Fruck'201 1). As at- 
mospheric attenuation is included in the shower reconstruction 



- devices responsible for telescope pointing 
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algorithm, air temperature and barometric air pressure fluctua- 
tions have an impact on the flux and the energy estimates of pri- 
mary gamma-rays. 

2.4. Barometric air pressure 

Generally, the air pressure analysis is performed by the astro- 
nomic community in order to see if a certain site is dominated by 
high pressure which would imply prevailing stable good weather 
(Lombardi et al. 2007 ). 

As one moves from the sea level to higher altitudes, the air 
pressure gets lesser and lesser. The theoretically expected baro- 
metric air pressure (Ptheo) at a given altitude can be calculated 
according to the US Standard Atmosphere model d 19761 1 with 
the following formula: 



P^Ph 



n 



Tb + Lb-{h- ht) 



(1) 



where go = 9,80665 m s" represents the sea-level value of the 
acceleration of gravity, Mq = 0.0289644 kg mol"' is the mo- 
lar mass of Earth's air, and R* = 8.31447 J K"' mol"' is the 
Universal gas constant. For altitudes lower than 11 km (tro- 
posphere), subscript b value is zero and it refers to sea-level. 
Then the geopotential height hb = 0, and the standard tempera- 
ture lapse rate Lb = -0.0065 K m"' . If we replace all these con- 
stants in Equation ([TJ, we obtain: 



F = Pol 1 -0.0065 — 



5.256 



(2) 



where Pq = 1013.25 hPa is the sea level standard atmospheric 
pressure, and To = 288.15 K (15 °C) is the sea level stan- 
dard temperature. For each investigated location, we have used 
Equation (|2]i to compute Ptheo and the results are included in 
Table [1] together with the geographical positions and the alti- 
tudes. The mean pressure value for a certain location, deter- 
mined from statistical analysis, is compared to the theoretically 
expected air pressure (Ptheo) for that altitude. The frequency with 
which it happens for the mean pressure to be higher than Ptheo is 
an indirect indication of high atmospheric stability (Jabiri et al. 
12000 ). On the other hand, variations of barometric air pressure 
in a short time scale (few hours) can induce weather instabilities 
with a serious impact on the operation of a telescope (Lombardi 
et al. 120091 ) ■ As the barometric air pressure may anticipate the 
temperature behavior by 2-3 hours (Lombardi et al. 120071 ). the 
information carried by the data can be used to better implement 
the thermalization of the telescope and of the instruments. In 
what concerns the Cherenkov telescopes, the barometric air pres- 
sure fluctuations impact on the atmospheric attenuation and im- 
plicitly on the shower reconstruction as discussed at the end of 
Subsection ' 



2.5. Sky cloudiness 

The clouds positioned directly overhead or not, can modify the 
response of a Cherenkov telescope. They can block the light 
coming from one direction and reflect the light coming from an- 
other This happens for the light generated within the extensive 
air showers, and also for the light originated from other sources 
(e.g. undesired city lights). This way either some useful light is 
lost or the sky may look brighter than it should be (Valcarcel 
120081 ). 



When clouds move into the field of view of a telescope 
they absorb Cherenkov radiation from extensive air showers and 
drops in the trigger rates result, rendering any background anal- 
ysis or absolute flux calculations very difficult (Cogan 2006)). 
A trigger rate lower than expected should always be accompa- 
nied by a lower number of recognized stars provided by the star- 
guider, in order to be sure that a cloudy atmosphere is the cause 
(Mazin 2007 ). In case only the trigger rate decreases, but the 
number of the recognized stars does not, this rather points to 
hardware problems of the telescope. 

Over the years, the standard procedure to evaluate the sky 
cloudiness has implied human observers going out, next to the 
telescope, and judging by themselves the condition of the sky. In 
the case of Whipple and VERITAS experiments, the observers 
gave a letter grade rating for the sky conditions at the moment 
of observations (Valcarcel l2008l Smith l2007l ). Even though, the 
procedure can be regarded as arbitrary because clouds are hard 
to see at night and it is considered inconsistent and subjective 
due to its human grading, it may still be regarded as a useful first 
criterion forjudging the validity of data (Smith 2007). 

An improvement of this situation was realized when a far- 
infrared (FIR) pyrometer connected to a telescope started to be 
used. The device was employed to detect the high thin wispy 
clouds that were difficult to see at night, by continuously mea- 
suring the temperature of the night sky (^elik 2008). This type 
of cloud could have on data, subtle effects which were poorly 
understood, and could have detrimental results on the analysis 
of source flux. 

Other tools for monitoring sky conditions included optical 
cameras (Kenny 120071 ) and weather satellite maps (Valcarcel 
20081 ). 

3. Statistical characteristics of the data 

For the present study we have sampled data from the archive 
of the National Meteorological Administration (NMA) of 
Romania. The meteorological data were collected in the period 
January 1, 2000 to December 31, 2009 by the NMA weather 
stations operated in the locations of Baisoara, Rosia Montana, 
Semenic and Ceahlau. The geographical positions and the alti- 
tudes of these locations are presented in Table[T]and more details 
can be found in (Radu et al. 120101 ). The weather stations were 
located in areas where reliable measurements and observations 
could be performed in accordance with the requirements of the 
World Meteorological Organization (WMO) (120081) . During the 
considered ten years period (2000 - 2009) no station was relo- 
cated. 

For comparison purposes, we have included in the present 
study, meteorological data recorded by weather stations de- 
ployed at the Nordic Optical Telescope (NOT) site. NOT is lo- 
cated at the Observatorio del Roque de Los Muchachos (ORM), 
in a region considered a reference for very good astronomical 
conditions - the Canary Islands, Spain. At ORM, NOT is po- 
sitioned in the proximity of MAGIC. The NOT meteorological 
data have been kindly made available to the public on the web 
page http://www.not.iac.es/weather/ 

In order to check for regional variations and long term trends 
we have analyzed meteorological data on relative humidity, dew 
point temperature, air temperature, wind speed, barometric air 
pressure, and sky cloudiness. The sky cloudiness parameter has 
been investigated only for the Romanian locations as there is no 
data acquired following a similar procedure in the case of NOT. 

The data records were available in the NMA database at reg- 
ular intervals of one hour, similar to the case presented in (Jabiri 
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Table 1. The geographical positions, the altitudes and the theoretically expected barometric air pressure (Ptheo) for the analyzed 
sites. 



Site 


Latitude 


Longitude 


Altitude [m] 


Ptheo [hPa] 


Rosia Montana 


46° 19' 03" N 


23° 08' 21" E 


1198 


877.5 


Baisoara 


46° 32' 08" N 


23° 18' 37" E 


1357 


860.6 


Semenic 


45° 10' 53" N 


22° 03' 21" E 


1432 


852.7 


Ceahlau 


46° 58' 39" N 


25° 57' 00" E 


1897 


805.1 


NOT 


28° 45' 26" N 


17° 53' 06" W 


2382 


758.0 



et al. I2000I I. As the NOT data were available at a higher fre- 
quency, for a proper comparison, we have first computed hourly 
averages. 

Nighttime (20:00 - 06:00), daytime (06:00 - 20:00) and en- 
tire day (00:00 - 24:00) intervals have been defined and all the 
data series are to be considered in local time. Monthly and an- 
nual samples have been extracted from these series of raw data 
according to a procedure to be described in the following. 

For a particular meteorological parameter (e.g. relative hu- 
midity), location (e.g. Baisoara), month (e.g. January) and time 
interval (e.g. nighttime) we have extracted data from the NMA 
database hour-by-hour, over the ten years period. The relevant 
statistics (median, standard deviation, skewness) have been com- 
puted for the obtained sample distribution. Subsequently, we 
have repeated the procedure for each of the 12 months and the 
obtained statistics have been used to produce a monthly distri- 
bution of medians. The procedure employed for the generation 
of the annual distributions of medians is similar, but instead con- 
sidering a certain month, data from a particular year (e.g. 2000) 
have been used. For the case of NOT, we have followed the same 
procedure, but hourly averages have been computed first. 

In addition to monthly and annual distributions of medians, 
we have also derived monthly and annual distributions of per- 
centages. This time, for a particular parameter, location, month 
and time interval, we have computed the percentage of data 
records whose values exceed certain limits (e.g. relative humid- 
ity > 90%). The procedure has been repeated for each month 
and the 12 resulted values have provided a monthly distribu- 
tion of percentages. The annual distributions of percentages have 
been obtained in a similar way. 

Within the univariate analysis performed, two statistical 
characteristics of each sample have been looked into: central ten- 
dency and dispersion. 

As measures of central tendency, arithmetic mean and me- 
dian are usually reported. Arithmetic mean is a useful statistic 
only if the investigated distribution is symmetrical (not too much 
skew), there are no outliers, and the analyzed data are measured 
at interval or ratio level. Median (the middle score in a set of 
scores) is useful when data are not symmetrically distributed or 
they are measured at an ordinal level. We have always reported 
the median statistic for the total sky cloudiness data which have 
been measured at ordinal level. In what concerns the other mete- 
orological parameters we have had to decide what the best statis- 
tic would be as a measure of central tendency. For each sample 
the amount of skewness has been computed using the formula: 



Vn(n — 1) m3 



n — 2 



(3) 



where m-}, - ^ TIi=i (.^i ~ ^ is the sample third central moment, 
"^2 = \ YIl=\ i^i ~ is the sample variance and x is the sam- 
ple arithmetic mean (Joanes & Gill 1998). However, the sample 
skewness tends to differ from the population skewness by a sys- 
tematic amount that depends on the size of the sample. In order 



to be able to tell, from the samples that we have, if the popu- 
lations are likely skewed (positively/negatively) or, if they are 
rather symmetrical, a test statistic is needed: 



(4) 



SES 



where SES is the standard error of skewness given by 



SES = 



6n (n - 1) 



i(«-2)(« + l)(«-i-3) 



(5) 



A population is considered to be very likely skewed if |Z| > 2 
(Cramer 1997). Otherwise no conclusion can be reached about 
the skewness of that population. It might be symmetric or 
skewed in either direction. 

We have computed this test statistic for all the samples used 
in the analysis. For most of the cases it has turned out that the 
populations have been very likely skewed. Therefore we have 
chosen to use as a measure of central tendency for all the dis- 
tributions involved in our analysis, the median statistic. This op- 
tion has no negative impact on those cases where the test statistic 
given by Equation (HJi, did not provide a clear conclusion. If the 
skewness is small, the difference between the median and the 
arithmetic mean is also small. For a symmetric distribution, the 
arithmetic mean and the median are equal. 

As an indicator for the dispersion of data, we have used the 
standard deviation, defined as the square root of the sample vari- 
ance m2. 

Following these considerations, the relevant statistics for all 
the data samples used in our analysis have been the median and 
the standard deviation. Whenever a plot in the present paper in- 
cludes error bars, the points represent the median values of the 
data samples for those specific months or years, while the error 
bars illustrate the standard deviation of the samples. As the er- 
ror bars frequently overlap, in order to have a clearer image of 
variability in the data, we have also included bottom plots of the 
standard deviation values. 

4. Astroclimatological results 

4.1. Relative humidity 

The relative humidity (RH) has been measured at the Romanian 
locations to accuracy of +2% by sensors placed 2 m above the 
soil surface. At NOT, the relative humidity sensors placed 2 m 
above the ground provided data to accuracy better than ±2% 
(Lombardi et al. 2007). 

In Figure[T] we present the monthly and the annual nighttime 
distributions of medians for relative humidity at the assessed 
sites. The bottom plots include the values of the standard de- 
viations for each month, and respectively for each year sample. 

The medians derived for the Romanian locations are clus- 
tered together and have values more than 40% higher in compar- 
ison to NOT. The driest Romanian site is Baisoara while Ceahlau 
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Fig. 1. The monthly (top-left) and the annual (top-right) nighttime distributions of medians for relative humidity. The bottom plots 
(left and right) show the values of standard deviations. 



is an extreme, unfavorable case displaying relative humidity me- 
dian values of 100% around the year (Figure [T] top-left and 
right). 

A seasonal dependence can be observed for the NOT data 
(Figure[T] top-left) and the lowest RH levels are recorded in July 
and August. A much fainter seasonal variation can be also ob- 
served for Rosia Montana and for Semenic. Baisoara exhibits 
an almost flat trend over the year, but from January to March, 
increased levels of relative humidity can be observed. 

The spread of data recorded at NOT exhibits a clear sea- 
sonal behavior (Figure [T] bottom-left) and it is larger than that 
of Romanian data (Figure [1] bottom-left and right). The mini- 
mum spread of NOT data is observed in July and August when 
the RH levels also reach a minimum (Figure [T] top-left). The 
spread of Romanian data shows a less visible seasonal variation 
with a minimum also observed during summer. 

The analysis of the annual distributions of medians shows 
that all sites exhibit flat trends over the years with fluctuations 
within the statistical errors (Figure [1] top-right). For all the 
Romanian locations but Ceahlau, 2000 appears as the driest year 
of the available data. 

In Figure|2] we display monthly and annual nighttime distri- 
butions of percentages for relative humidity data records whose 
values are larger than 90%. Considering this value as a safety 
operational limit, the points in this figure indicate the percent- 
age of time that a telescope built at the considered locations 
would have been unusable because of high levels of RH. The 
larger are these percentages the worse are the observational con- 
ditions. A seasonal variation is visible for all locations (Figure|2] 
left) and the summer months appear to be best suited for obser- 
vations everywhere but at Ceahlau where a small improvement 
can be observed only from October to December. The small- 
est computed values belong to NOT and they confirm that its 
site is excellent for astronomical observations. Among the in- 
vestigated Romanian locations, Baisoara offers the closest con- 
ditions to those of NOT, especially during October, November 
and December. 



In what concerns the annual nighttime distributions of per- 
centages for RH (Figure |2] right), Baisoara shows close values 
to NOT, but an increasing trend can be observed since 2006. This 
raises the question if we are in the presence of a change in the 
local microclimate or it is a typical oscillating behavior of the 
relative humidity over the years. More data are necessary to set- 
tle the issue. 

Figure[3]presents the monthly and the annual nighttime, day- 
time and entire day distributions of medians for relative humid- 
ity at Baisoara, Rosia Montana, Semenic and NOT. The median 
values of data recorded at Ceahlau are constantly equal to 100% 
and not included in this figure. The nighttime median values are 
higher than the daytime ones (Figure |3] top and bottom) and 
night-day variations of ~ 20% can be observed especially dur- 
ing summer (Figure[3] top) at the three Romanian locations. The 
lower levels of RH during daytime can be exploited for main- 
tenance operations. At the NOT site, the night-day variations 
are very small and sometimes the nighttime median values are 
smaller than those of daytime data. This is especially favorable 
for the astronomical observations as they take place at night. 



4.2. Dew point temperature 

In order to evaluate the risk of condensation, we have analyzed 
the monthly and the annual nighttime distributions of percent- 
ages for the difference between dew point temperature and air 
temperature data records, whose values are smaller than an up- 
per limit. The larger are the values of these percentages for a cer- 
tain site the higher is the risk of condensation expected for that 
particular site. Two values have been considered in our study for 
the upper limit: 1 °C and 5 °C. For the analyzed meteorologi- 
cal data of all five locations, the cases where this difference is 
smaller than 1 °C also imply relative humidity values larger than 
90%. As this situation has been analyzed in Subsection ^4.11 we 
present here only the case when the upper limit is set to 5 °C 
(Figure HJ. This can be regarded as an extension of the analysis 
performed in Subsection 34.1l to values of RH lower than 90%. 
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Fig. 2. The monthly (left) and the annual (right) nighttime distributions of percentages for relative humidity data records whose 
values are larger than 90%. 
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Fig. 3. The monthly (left) and the annual (right) nighttime, daytime and entire day distributions of medians for relative humidity. 
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The location characterized by the lowest risk of condensa- 
tion is that of NOT. The statistical distribution of data taken at 
its site shows a clear seasonal dependence (Figure [H /e/f) with 
a minimum reached in July and August. A seasonal dependence 
is also visible for Rosia Montana and Semenic with lower val- 
ues reached during summer. Baisoara exhibits the lowest risk of 
condensation during April and May, while high values have been 
observed from January to March. Ceahlau is again the worst case 
even if for the interval October - December a decrease for the 
risk of condensation can be noticed. 

The analysis of the annual distributions (Figure |4] right) 
shows a dry climate at NOT and confirms Baisoara as the dri- 
est Romanian location. The Romanian sites exhibit an oscillat- 
ing behavior over the years and 2000 appears as the least humid 
one. 

4.3. Wind speed 

At the Romanian locations, wind speed has been measured by 
sensors placed 10 m above the soil surface, to accuracy that de- 
grades as speed increases. The NOT wind sensor took data to 
accuracy better than 2% (Lombardi et al. 120071) . 

In Figure |5] we show the monthly and the annual nighttime 
distributions of medians for wind speed at the assessed sites. A 
clear dependence on altitude (see also Table\I} can be noticed for 
the median values and for the spread of data. The higher is the 
altitude of a location, the stronger are the wind speeds and the 
larger are the variations of data. Baisoara is a favorable exception 
as its site has the lowest wind speeds even if it is not located at 
the lowest altitude. This effect is probably due to the local orog- 
raphy. The unfavorable exception is Ceahlau which has similar 
wind conditions to NOT and a larger spread in data, even if it is 
situated lower than that. 

Faint seasonal variations can be observed for the distribu- 
tions derived at NOT, Ceahlau, and Semenic (Figure |5] top-left) 
and the lowest wind speeds are measured during summer A sea- 
sonal variation is also visible for the spread of all Romanian data 
(Figure|5] bottom-left). This dependence is clearly observable at 
Ceahlau, but it is weaker for the other Romanian locations. 

The analysis of the annual distributions (Figure|5] top-right), 
shows that all sites exhibit almost flat trends over the years with 
fluctuations within the statistical errors. 

Considering 15 m i ' as the wind speed Umit beyond which 
a telescope should be turned off and brought to the parking posi- 
tion ( ^2.2l i. we have plotted in Figure |6]the monthly and the an- 
nual nighttime distributions of percentages for wind speed data 
records larger than this limit. The larger are these values for a 
certain site the higher are the percentages of time when a tele- 
scope operated there should be turned off because of high winds. 
All the Romanian locations, but Ceahlau show values lower than 
the NOT site and almost flat trends. The distribution of Ceahlau 
data shows a seasonal dependence (Figure|6] left) with the lowest 
values from April to September when its conditions are closer to 
those at NOT. 

A good picture of the wind speed distributions at the ana- 
lyzed locations can be obtained from Table |2] It reports for the 
period April - October, over the ten years range, the percentages 
of nighttime wind speed data records whose values fit into fixed 
intervals. Baisoara has a predominance of wind speeds smaller 
than 3 m (71.5%) whereas the other locations behave worse. 
However, this parameter exceeds 10 m i"' only in 0.4% of cases 
at Rosia Montana and respectively 5.8% at Semenic, so these 
locations also preserve good wind speed conditions. 



Table 3. Nighttime duty-cycles computed at Baisoara (BA), 
Rosia Montana (RM), Semenic (SE), Ceahlau (CE) and NOT, 
as the percentage of nighttime data records with RH < 90% and 
wind speed < 15 m ' , from April to October (2000 - 2009). 



BA 


RM 


SE 


CE 


NOT 


Duty-cycle [%] 78.9 


68.6 


56.4 


16.5 


87.9 



In case 15 m s is considered as a safety operational 
limit, a telescope installed either at Baisoara, Rosia Montana, 
or Semenic will practically never have to be turned off because 
of strong winds. 

If in addition to this wind speed limit we also consider that 
the relative humidity values should not exceed 90%, a duty-cycle 
can be computed for each analyzed location. We have calculated 
this parameter only with nighttime data recorded from April to 
October over the ten years period and the results are presented 
in Table |3] 

The best result for a Romanian location has been obtained 
for Baisoara that would offer a duty-cycle of ~ 79%. This value 
is quite close to the one obtained for the reference site of NOT 
(~ 88%) and it recommends Baisoara as a good potential site for 
the establishment of a Cherenkov telescope in Romania. 

4.4. Air temperature 

The air temperature at the Romanian locations was measured by 
sensors placed 2 m above the soil surface to accuracy better than 
+0.2 °C. The same type of measurements was performed at NOT 
to accuracy of ±0.1 °C. 

In Figure |7] we display the monthly and the annual night- 
time distributions of medians for air temperature. A clear sea- 
sonal variation is noticeable for all investigated sites, but the 
difference between winter and summer months is higher for the 
Romanian locations (Figure|7] top-left). Ceahlau is characterized 
by the lowest air temperatures around the year, while the other 
Romanian locations exhibit similar median values. From April 
to October the air temperature conditions at Baisoara, Rosia 
Montana and Semenic are very similar to those at NOT. 

A larger spread characterizes air temperature data at the 
Romanian locations in comparison to NOT (Figure |7] bottom- 
left and right) and a clear seasonal dependence can be observed 
from Figme^(bottom-left). Less variability in the Romanian air 
temperature data occurs from April to September when its values 
tend to get closer to those of the NOT site. 

The analysis of the annual distributions (Figure |7] top- right) 
shows for all sites almost flat trends with fluctuations within the 
statistical errors. 

Figure |7] (top-left and right) also depicts an altitude depen- 
dence (see also Table \J} for the median values of air temper- 
ature distributions. Ceahlau, the highest altitude of the investi- 
gated Romanian locations is also the coldest one. As the alti- 
tude decreases, the temperature values increase. Rosia Montana 
and Baisoara exhibit the highest temperatures. Even if located at 
the highest altitude, NOT is also characterized by high temper- 
atures. This effect can be explained by the much lower latitude 
where NOT is positioned in comparison to the Romanian loca- 
tions (Table [TJ. The increase in altitude is compensated by the 
decrease in latitude and the median temperatures increase. 

Special technical procedures can be employed to overcome 
the negative effect of freezing on the proper functioning of tele- 
scopes. However, an ideal site to host a telescope would very 
rarely have air temperatures below °C. Figure [8] shows the 
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Fig. 4. The monthly (left) and the annual (right) nighttime distributions of percentages for the difference between dew point temper- 
ature and air temperature data records, whose values are smaller than 5 °C. 
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Fig. 5. The monthly (top-left) and the annual (top-right) nighttime distributions of medians for wind speed. The bottom plots (left 
and right) show the values of standard deviations. 



monthly and the annual nighttime distributions of percentages 
for air temperatures data records whose values are below °C. 
The larger are these values for a certain site the higher are the 
periods of time when a telescope operated there would be af- 
fected by freezing. From the analysis of these distributions re- 
sults that in Romania, the best period for observations is May - 
September (Figure [8] left) when data are quite similar to those 
at NOT and there is practically no freezing. The annual distribu- 
tions (Figure[8] right) show larger year-to-year fluctuations at the 
Romanian sites in comparison to NOT (exception for 2000 and 
2001). Baisoara and Rosia Montana experience the least freezing 
among the Romanian sites. 



The median values of the monthly and the annual nighttime, 
daytime, and entire day distributions of air temperature are pre- 
sented in Figure|9] For all locations, the largest daytime to night- 
time variations can be observed during summer (Figure |9] top) 
and larger differences can be noticed at the Romanian locations 
compared to NOT. The annual distributions (Figure |9] bottom) 
show similar levels for the daytime temperatures at Baisoara, 
Rosia Montana and NOT, while the nighttime temperatures are 
lower in Romania. 
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Fig. 6. The monthly (left) and the annual {right) nighttime distributions of percentages for wind speed data records whose values 
are larger than 15 m i"' . 

Table 2. Nighttime wind speed statistics at Baisoara (BA), Rosia Montana (RM), Semenic (SE), Ceahlau (CE) and NOT, collected 
in a ten years time period (2000 - 2009), from April to October. 



tV,,p range [m BA [%] RM [%] SE [%] CE [%] NOT [%] 

< Wsp < 3 
3 < W,p < 10 

^ sp 

> 15 



71.5 


52.6 


28.9 


23.5 


16.5 


28.0 


47.0 


65.0 


51.1 


65.3 


0.5 


0.4 


5.8 


17.9 


15.4 


0.0 


0.0 


0.3 


7.5 


2.8 



4.5. Barometric air pressure 

The barometric air pressure was measured at the Romanian loca- 
tions 1.20 m above the soil surface to accuracy of +0.1 hPa. The 
NOT air pressure sensor was placed 2 m above the ground and it 
provided data to accuracy of ±0.1 hPa (Lombardi et al. 12007b . 

Figure [To] shows the monthly and the annual nighttime dis- 
tributions of medians for barometric air pressure at the studied 
sites. A clear inverse correlation with altitude can be observed 
in the top-left and the top-right sections of the figure {see also 
Table\J}. The highest altitude site (NOT) exhibits the lowest air 
pressures, whereas the highest pressures have been recorded at 
the lowest altitude site (Rosia Montana). 

A very faint seasonal variation can be observed from 
Figure [TO] {top-left) at all Romanian locations with higher me- 
dian values of barometric ak pressure recorded during summer. 
The analysis of the annual distributions (Figure (TO] top-right) 
shows an almost constant flat trend over the years for all loca- 
tions. 

The spread of data recorded at the Romanian sites shows a 
more pronounced seasonal dependence and it is larger in com- 
parison to NOT (Figure [TOl bottom-left). The lowest spread of 
data is observed during summer when the weather conditions 
are more stable and the air pressure values fluctuate less. The 
lower spread of NOT air pressure data is also an indicator of a 
more stable weather. 

When the local barometric air pressure is lower than the the- 
oretical value (Ptheo) at the same altitude, the weather conditions 
may become unstable ( ^2. 4b . In order to compare the stability 
of weather conditions, we have computed for the investigated 



sites the percentages of local barometric air pressure data records 
whose values are lower than those theoretically expected. The 
lower are these percentages for a certain site, the better are the 
observational conditions offered by that particular site. The re- 
sults are presented in Figure [TT] 

The percentage values at NOT are roughly equal to 0%. 
Almost always the values of local barometric air pressure are 
higher than those theoretically expected. The site is dominated 
by high pressure values and remarkable stable weather condi- 
tions characterize it. The situation is different at the Romanian 
sites where the values of the local barometric air pressure often 
fall below those theoretically expected. Figure [TT](/e/f) shows a 
very clear seasonal dependence with the most favorable weather 
conditions for astronomical observations from May to October. 

The analysis of the annual distributions (Figure [TT] right) 
shows that for all Romanian sites these percentages fluctuate ex- 
tensively. Large differences can be noticed for the years 2000, 
2001, 2002, 2004, and 2007, while smaller ones have appeared 
in 2003, 2005, 2006, 2008, and 2009. 

From the distributions depicted in Figure [TT] Semenic ap- 
pears to have the most stable weather of the Romanian investi- 
gated sites. 

4.6. Sky cloudiness 

For our analysis we have used a parameter stored in the database 
of the NMA and called "total sky cloudiness". This is a param- 
eter characteristic for the sky condition and it represents the de- 
gree to which the sky is covered by any type of clouds at the 
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Fig. 7. The monthly (top-left) and the annual (top-right) nighttime distributions of medians for air temperature. The bottom plots 
(left and right) show the values of standard deviations. 
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Fig. 8. The monthly (left) and the annual (right) nighttime distributions of percentages for air temperatures data records whose 
values are below °C. 



moment of observation. For the Romanian locations and the pe- 
riod of time considered (2000 - 2009) this parameter was vi- 
sually determined by specially trained human observers. They 
went out on meteorological platforms and from a spot that al- 
lows maximum visibility of the celestial vault, evaluated the to- 
tal sky cloudiness on a scale that stretched from to 9. "Zero" 
meant a "perfectly clear sky". "Eight" stood for a "completely 
covered sky" and "nine" was a "sky invisible because of fog". 

Figure [T2]exposes the monthly and the annual nighttime dis- 
tributions of medians for total sky cloudiness at the Romanian 
sites. The monthly nighttime distributions of total sky cloudiness 
(Figure [12] top-left) show a seasonal variation with the lowest 
cloudiness during summer Rosia Montana and Baisoara offer 



the lowest levels of cloudiness and August appears to be the best 
month for observations. The analysis of the annual distributions 
(Figure[T2] top-right) confirms the quaUty of these two locations. 

Baisoara and Rosia Montana show a seasonal dependence 
for its spread of data with the lowest values obtained during sum- 
mer (Figure [12] bottom-left). The standard deviations included 
in the plot of the annual distributions (Figure [12] bottom-right) 
show the lowest spread for the data collected at Baisoara. 

The levels of total sky cloudiness data recorded in the NMA 
database do not reflect the thickness of the cloud layers and 
implicitly the quality of the atmospheric transmission. A level 
eight sky is a completely covered sky, but the layer of clouds 
can be a thin one and the VHE gamma-ray observations might 
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Fig. 9. The monthly (top) and the annual (bottom) nighttime, daytime and entire day distributions of medians for air temperature. 
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Fig. 11. The monthly (left) and the annual (right) nighttime distributions of percentages for local barometric air pressure data records 
whose values are lower than the theoretically expected barometric air pressure values given in Table [1] 
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Fig. 12. The monthly (top-left) and the annual (top-right) nighttime distributions of medians for total sky cloudiness. The bottom 
plots (left and right) show the values of standard deviations. 



still be possible. On the other hand, under a level five sky (partly 
covered), thick clouds may frequently move into the telescope's 
field of view and it may seriously impede on taking high quality 
data. 

Under these circumstances we have computed the percent- 
ages of total sky cloudiness data records whose values are 
strictly smaller than 2. A sky of level or 1 is perfectly clear 
and it will certainly allow high quality observations of very high 
energy gamma-rays. The higher are these values for a certain site 
the better are the observational conditions offered there. The re- 
sults are presented in Figure [13] (left and right) and they show 
that Baisoara and Rosia Montana exhibit the best observational 



conditions, while August is the month when sky is most fre- 
quently clear. 

In Figure [14] we have plotted the monthly and the annual 
nighttime and daytime distributions of medians for total sky 
cloudiness. At Baisoara, the total sky cloudiness is lower dur- 
ing nighttime than during daytime, around the year This is a 
particularly favorable eff'ect for the astronomical observations as 
they take place at night. For the other locations (Figure [14] top), 
the effect is observable only on certain periods (e.g. March - 
September at Ceahlau). The annual distributions (Figure[T4] bot- 
tom) show nighttime lower than daytime values of total sky 
cloudiness at Baisoara and Rosia Montana over the analyzed ten 
years period. 
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Fig. 13. The monthly (left) and the annual {right) nighttime distributions of percentages for total sky cloudiness data records whose 
values are strictly smaller than 2. 



(I) 

« 
C 

u 



in 

ID 

« 

3 

o 

CO 





— Nighttime values 
Day time values 
















\-\ ,-" ' / 






^^^^^^^^^^^^^^^^^^^^^^^^ 






BAISOARA 


ROSIA MONTANA 

. 1 . 1 . 1 . 1 . 1 . 1 


SEMENIC y 

. 1 . 1 . 1 . 1 . 1 . 1 


CEAHLAU 

1 . 1 . 1 . 1 . 1 . 1 


2 4 6 8 10 12 


2 4 6 8 10 12 
Month 


2 4 6 8 10 12 


2 4 6 8 10 12 




— Niglit time values 








— Day time values 








/■"'\ / \ / 














BAISOARA 


ROSIA MONTANA 


SEMENIC 


CEAHLAU 











20002002200420062008 2000 2002 2004 2006 2008 2000 2002 2004 2006 2008 2000 2002 2004 2006 2008 

Year 

Fig. 14. The monthly (top) and the annual (bottom) nighttime and daytime distributions of medians for total sky cloudiness. 
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5. Conclusions 

A site assessment campaign has been initiated for the identifi- 
cation of the best location where to install and operate a small, 
imaging atmospheric Cherenkov telescope in Romania. 

The results of an astroclimatological study based on the 
statistical analysis of relative humidity, dew point temperature, 
air temperature, wind speed, barometric air pressure, and sky 
cloudiness local data have been presented in this paper. 

Data over a ten years period (2000 - 2009) have been 
sampled from the archive of the National Meteorological 
Administration (NMA) of Romania, for four candidate locations 
(Baisoara, Rosia Montana, Semenic, Ceahlau). Meteorological 
data recorded at the Nordic Optical Telescope site (Canary 
Islands) have been used as a reference for comparison purposes. 

Prior to presenting the results of the analysis, the influence 
of atmospheric conditions on the operation of Cherenkov tele- 
scopes has been reviewed and the statistical characteristics of 
the available meteorological data have been studied in order to 
choose the relevant statistics of distributions. 

The analysis of relative humidity data has revealed, as ex- 
pected that the weather in Romania is more humid than at NOT. 
However, among the Romanian investigated sites, the microcli- 
mate of Baisoara appears to be the driest one and it offers the 
closest conditions to those reported for NOT. During summer 
the relative humidity nighttime data exhibit the least spread and 
they are clustered around low level central values at all inves- 
tigated locations, but Ceahlau. Night-day variations of up to ~ 
20% have been observed between the medians of the statistical 
distributions at the Romanian sites, while at NOT these varia- 
tions are very small. 

In order to evaluate the risk of condensation, we have an- 
alyzed the nighttime distributions of percentages for the differ- 
ence between dew point temperature and air temperature data 
records, whose values are smaller than 5 °C. The lowest risk 
of condensation is characteristic to NOT, while in Romania, 
Baisoara offers the best conditions. 

The wind speed analysis has shown that the higher is the 
altitude of a site, the stronger are the winds and the larger is 
the spread of data. Baisoara is an exception being characterized 
by the lowest winds, even if it is not located at the lowest alti- 
tude. Speeds smaller than 3 m s ' have been measured there in 
71.5% of the time. Wind speed conditions are generally better in 
Romania than at NOT especially due to the lower altitudes. 

A duty-cycle has been computed considering that a telescope 
should be turned off in case wind speed exceeds 15 m s"' or 
the relative humidity is larger than 90%. For samples including 
nighttime data recorded from April to October over the ten-years 
period, the best result of the Romanian locations has been ob- 
tained for Baisoara (~ 79%). This value is quite close to that 
obtained for the reference site of NOT (~ 88%). 

Very similar air temperature conditions have been noticed at 
Baisoara, Rosia Montana and Semenic in comparison to NOT 
for the period April to October and practically no freezing oc- 
curs at these locations from May to September An inverse cor- 
relation has been observed between air temperature and altitude. 
Rosia Montana and Baisoara exhibit the highest temperatures, 
while Ceahlau is the coldest one. Even if located at a higher al- 
titude than Ceahlau, NOT exhibits high temperatures. This ef- 
fect can be explained by the much lower latitude where NOT 
is positioned in comparison to the Romanian locations. Slightly 
larger night-day variations of air temperature have been detected 
in Romania during summer compared to NOT. 



The analysis of the barometric air pressure data has pointed 
out that the NOT site is dominated by high pressure values 
and remarkably stable weather conditions characterize it. At the 
Romanian sites the weather conditions are more stable from May 
to October and Semenic has been designated as the location with 
the most stable weather. The inverse correlation between the 
barometric air pressure and the altitude has been observed for 
all locations. 

Data of total sky cloudiness have been available only for the 
Romanian locations and our study has indicated that Baisoara 
and Rosia Montana offer the largest percentages of time when 
the sky is perfectly clear The analysis of the night-day varia- 
tions has allowed us to observe that at Baisoara the nighttime 
cloudiness is lower on average than at daytime. This effect is 
particularly favorable for the astronomical observations. 

The results of the present study have made us select Baisoara 
as the best investigated site for the establishment of a Cherenkov 
telescope in Romania. Good meteorological conditions have 
been also determined at Rosia Montana and as these two sites 
are both located in the Apuseni Mountains, we consider this area 
as the optimal place for the location of the telescope. A suitable 
period, around the year, for operating a Cherenkov telescope in 
Romania, appears to be from April to October, but taking data 
is not restricted to this time frame only. However, beyond this 
period the number of usable nights is lower 
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